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Abstract 

In the present work, we study the classical behavior of an electric dipole in presence of an 
external uniform magnetic field. We derive equations and constants of motion from the Lagrangian 
formulation. We obtain an infinitely periodic effective potential that describes a rotational motion. 
The problem is not directly separable in relative and center of mass variables; even though, we 
are able to write the energy of the system as a function of an only term, the relative variable. We 
define another constant of motion, which couples the relative with the center of mass variables. 
We describe conditions for bound states of the dipole. In addition, we discuss the problem in 
the approximation of small oscillations. Finally, we explore the existence of a possible family of 
trapped states in a region of the space where there are no classical turning points. 
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I. INTRODUCTION 



take the 



Nowadays, fundamental elements in scale of molecular machinery 
attention from several specialists, due to the development of nanotechnology. An alternating 
electric field has already used to explore electronic structures, however this process can be 
useful to exchange the orientation of molecules; therefore, it is possible to obtain some 

. Devices on the molecular level 
are obtained from the conversion of energy into controlled motion; regardless of this, it is 
difficult to repeat this process in a mechanical molecular motor. This process is common in 
biological systems Q]. By moment, it is expected to find physical principles of a motor in 
molecular scale through the dynamics of rotors in two dimensions. Such rotors are modeled 
as electric dipoles in electric or magnetic fields. 

It is well known, we know that the simplest many-body system is the two-body system and 
we expect to separate the two-body problem in two problems of a single body. The problem 
of two equal charges in an external magnetic field fulfills to this requirement completely. In 
a suitable formulation, two independent variables are considered to describe the motion of 
the center of mass and the relative motion, respectively. The classicaljsl and quantum |lo| 
physics of the problem are clearly described in details . On other hand, the problem of two 
particles with opposite charges into a magnetic field does not allow to make the separation 
of variables in independent equations, as it is possible to do it in the case presented above. 
Nevertheless, two constants of motion are obtained from a suitable Lagrangian formulation. 
It is possible to write an equation for the energy in terms of the relative variable and other 
constant of motion connecting the relative variable with the center of the mass coordinate. 
Cla.ica.a and ,ua,.u.H result, were previously .epo«ed. 

In addition, other related cases are discussed in precedent works 

• the dynamics of charges in oscillating electric and magnetic fields in Ref . [3, 0] , 

• the dynamics of a dipole in a magnetic field on Y direction, the motion of its center 
of mass is restricted to the Z direction, and its rotation to the XY plane, is discussed 
in Ref.Q- 

The main goal in this work is to describe the classical behavior of an electric dipole in an 
external magnetic field. With this purpose, we introduce a model for an electric dipole 
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as two particles of arbitrary masses, where separation between them is constant, and its 
charges are equal and opposite. That case is a two-body problem and its exact analysis 
requires a difficult treatment and the further useful understanding, which is obtained by 
approximations. 

The paper is organized as follows. Firstly, we present a general formulation for an electric 
dipole in magnetic field. We define the Lagrangian function and we derive constants and 
equations of motion. Secondly, we present the solution for a particular problem; this is, the 
planar motion of an electric dipole in a perpendicular magnetic field. We discuss bound and 
trapped states. Finally, we summarize the formulation and main results. 



II. LAGRANGIAN FORMULATION: SYMMETRIC GAUGE 

In the present model, internal coupling holds the two charges of the dipole together and 
the Coulomb interaction between the charges comes to be constant. Then, we consider a 
rigid dipole, two fixed charges by a massless rod, in the presence of a uniform magnetic field. 
One of particles carries charge +e, while the other — e. The magnetic field is derived from a 
vector potential A, as follows B = V x A. We relate to the particle 1(2) the position ri(r2), 
the velocity ri(r2) and the mass 1711(7712) . The Lagrangian formulation leads to the following 
expression: 

L{fi, r2] ri, f2) = -niir^ + -m2r2 A{ri) ■ fi + -A{f2) ■ r2 + —r^ (1) 

2 2c c K\r2 — Til 

where k is the dielectric constant of the medium in which the motion of the dipole occurs. 
We define the vector potential A through the symmetric gauge as follows 

A{n) = \bx f,, for 1=1, 2, (2) 

— * 

where B is the uniform magnetic field. Now, we take into account the following change of 
variables 

f=f2-fi, (3) 

- mifi + m2r2 

^ = : ) [V 

ITLi + m2 

where r is the relative position and R is the position of the center of mass. Now, if we 
replace Eq.(j2I)-Eq.(jH)into Eq.((H), we obtain the following function 

^ • 1 ^2 X -2 

L{R,f;R,r) = -MR + -fir + 



2 ' 2^ ' K\f\ 



e 

+ - 

c 



-BxR-f+-Bxf-R+— -^-Bxf-f\ (5) 

2 2 M 2 ^ ^ 



mass. 



where M = mi + m2 is the mass of the center of mass and /x = 71111712/ M is the reduced 
From now on, we consider mi = m2. In contrast to the case 
the present Lagrangian function is similar to particles with opposite charges because it is a 
nonseparable problem; the motion of the center of mass is coupled to the relative variable j^. 
The conjugate momentum for the center of mass is given by 



= ^ ' V ' ^ = MR + ^B xr, (6) 



Pr = ^ ' ^ = -i-J X -f (7) 



dL{R, f]R,r) J^^^ , e 
OR 

which depends on center of mass and relative motions. From deriving the equations of 
motion we have 

dL{R,f;R,f') _ e 
'dR ~ ~^c' 

If we integrate the Eq.(|7j), and it compares with the Eq.®, we obtain the first constant of 
motion 

MR + -B xr = C, (8) 

c 

and the first equation of motion is given by 

MR = --B X 'r. (9) 
c 

We obtain for the relative vector conjugate momentum 



dL{R,f;R,r} 1. e - - 

Pr = = yur + —B X R, (10) 

or 

thus, the force is given by 

u dL{R,r;R,r) e - 5 „ 

Pr = = — B X R- —^Cr. 11 

or c K\f\^ 

If we derive the Eq. pO|) in time, and we compare with Eq. pH) . we can obtain the second 
equation of motion 

lj,r=--BxR ^e;. (12) 

C K\f\ 

The second constant of motion is E, the energy of the system. The general way to calculate 
the energy from a knowledge L is 

E = E{R, f;R,r) = Pji- R + pr-f- L{R, r; R, f), (13) 



and if we take Eq.(jHI) and we replace it into Eq. ()13p we can explicitly obtain 



E = iMif + i,,^ - (14) 

Putting Eq.(jHl) in Eq. p4j) . we obtain an expression for the energy in terms of the relative 
variable only. 

1 „ 1 I „ 2 „2 

(15) 



li^ 1 \e ^ ^ ^ ^ e 



E = -nr \-Bxr — C 



K\f\ 



A particular case has been discussed in Ref.jl^, where the motion of the center of mass of 
the dipole is restricted to the Z direction and rotational degrees of freedom of the dipole to 
the XY plane and the magnetic field is confined to the Y direction. Hence, we recover the 
classical behavior of the system from basic equations of the present formulation. 



III. MOTION IN A PERPENDICULAR PLANE 



We will choose the most natural motion of electric charges in a magnetic field. Let us 
consider a confined motion to two dimensions in a perpendicular plane to the magnetic field. 
In this way, we restrict the degrees of freedom of the center of mass and the relative motion to 
the plane XY, the direction of the uniform magnetic field is chosen to be in the Z direction, 
this is, B = Bz. The separation between charges is constant, and |r| = a is the linear size 
of the dipole. Let us define variables and parameters in terms of the dimensionless units: 
p = f/a , ^ = R/a, p = r/tOcO-, ^ = R/uca, if Eq = ^Ma^cUc^, then e = E/Eq, Lpc = f^/ujc, 
9c = 9/uJc , where Uc = eB/Mc is the cyclotron frequency and a is the size of a dipole. 
We define Ec = e'^/Eona as the ratio between the Coulomb energy and magnetic energy. In 
summary, distances have been defined in terms of a and the time in term of 1/wc- 

From Eq.(jHl) and considering the dimensionless parameters, we define a dimensionless 
constant of motion as follows 

C,+zxp = zxT, (16) 

— * 

where T is a vector constant of motion and p represents a position in the perpendicular 
plane XY to the magnetic field. In addition, from Eq. ()14|) and considering the dimensionless 
parameters, we write the energy as a dimensionless equation 

e=\^n\\'p\'-ec. (17) 
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FIG. 1: An illustrative picture is shown for the angle ^. In addition, other parameters related to 
the geometry and the orientation of the dipole are depicted. 

Now, if we replace Eq. fll6p into Eq. p7|) . we write the energy as a function of a single variable; 
this is, 

e = \\'p\'+\p-f\'-e,. (18) 
The dimensionless effective potential is given by 

K// = |p-fp-e,. (19) 

By following, if the linear size of the dipole is the unit of the distance (i.e., |p| = 1), we can 
write 

Veff = 1 + f2 - 2T COS ^ - Ec, (20) 

where represents the angle between the line that joins the charges of the dipole and the 
vector constant T. Thus, the effective potential Ve// is a periodic function of Finally, we 
write an explicit expression for the dimensionless energy given by 

5= 1^2^1 + f2-2TcOS^-£e- (21) 

In addition, in the Figure (^, angle dimensionless relative variable and the velocity of 
the center of mass are depicted at fixed time. Several parameters related to the geometry 
and the orientation of the dipole are shown. 
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A. Bound states 



According to previous considerations, we know that the effective potential is a periodic 
function of the angle If the energy e of the system lives between two extreme values, 
V^Y = ((1 ~ '^)^ ~ ^c) and V^'j^ = ((1 + T)^ — Sc)-, the motion is bounded by two turning 
points. The value of total energy in such point e = Veff{p) gives the limit of the oscillation. 

Hence, from the Eq. ()21|) and by considering that p ■ 4/ = d'^/dr, we obtain by direct 
integration 



Ar(^) =r(^) -r(^o) 



1 



(22) 



2J^o[e + e^-l-T^ + 2T cos ^']i/2 ' 
where r = ujct, which we write in terms of an incomplete elliptic integral of the first kind as 



follows 



Ar(M/) = -F 



arcsm 



' sin\l//2 ' 
sin \l>£/2 



^ 2 



(23) 



where r = Uct 

The typical oscillation of the variable \1/ for bound states represents a bound rotation 
between two limiting angles, ±\I^£, the classical turning points. That limiting angles are 
given by 



,e + Sc 



arccos 



2T 



(24) 



From the Eq. ()18|l . the energy of the system determines the amplitude of the the oscillation. 
According to the Eq. ljlUp . if the dipole enters to the magnetic field without an angular 
relative velocity, then the energy only depends on the initial orientation of the dipole. From 
the Eq.((221) and the Eq . ()24|) the period of the motion can be given by 



Tp 



d^' 



(25) 



/2T J-^e [cos ^' - cos ■ 

Now, we use the definition of the complete elliptic integral of the first kind to write the 
period of the motion 

TP = -^K [sin ^1 . (26) 

n 

In addition, we can use the definition of the hypergeometric function[lj| to present the same, 
as follows. 



rp 



2^1 



1 1 



- , - , sm — 
2' 2' 2 



(27) 
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FIG. 2: The frequency is depicted as a function of energy e + Ec- The exact value for the frequency 
coincides with the value in the aproximation of small oscillations in the limit e —£c 

Using the approximation of small oscillations, i.e., for energies near to the point of the stable 
equilibrium V^T", the dimensionless frequency up is given by 



Clearly, this frequency depends on T, which corresponds to the initial orientation and the 
initial velocity of the center of mass from the Eq. ()ll|) . 

In the Figure El is depicted the relation between the dimensionless frequency and the 
dimensionless energy. In the approximation of small oscillations, we can get the effective 
potential as a quadratic function of the angle \1'. This is, the period of such oscillations is 
independent of its amplitude. The exact value for the frequency u foi e + Ec ^ approaches 
to the analytical value for the frequency ujp from small oscillations. 



(28) 



UJp = 



2 
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B. Trapped States 



The classical trapped states are defined when the mean value of the center of mass velocity 
is zero. From Eq. (jl6|) we obtain 

(f) = - r dr'i-z X p + z X f) = (29) 
r Jo 

where r is the dimensionless period. The present condition represents a closed motion. We 
expect to characterize a family of unbounded states, without classical turning points, which 
can satisfy the previous restriction. By following, we write by compounds, 

{(u) = -- r dr' COS 9 = (30) 

T Jo 

{f±) = - r dr'i- sin ^ + T) = 0, (31) 

r JO 

where 6 is the angle between ^ and T (see Figure P) and the symbols and are referred 
to the parallel and perpendicular parts of ^ respecting to the vector T. Furthermore, we 
take into account that T = TT, where T is a unit vector. We see that if T = 0, Eq. ljnnj) and 
Eq. (jHT|) are simply satisfied, the classical motion is trapped for any dimensionless energy 
e > 1 — we obtain from Eq.(|2ip the frequency of the trapped state that is given by 
\[' = 2y/6 + Ec — 1- In addition, we can emphasize that if T 7^ 0, Eq. (jHn|) is always satisfied; 
however, the Eq. (jHT|) imposes that — 1 < T < 1, while if |T| > 1 the trapping should not be 
possible. 

IV. SUMMARY 

We have proposed to introduce a description of the classical dynamics of rotors such 
as electric dipoles in an external magnetic field. Our treatment does not consider radiation 
effects; however, the solution is nontrivial; therefore, we believe the present discussion would 
be pedagogically useful as a special topic of a standard course of classical mechanics. 

In the present application, we have considered the motion in a perpendicular plane to a 
uniform magnetic field. The motion axis coincides with the direction of the magnetic field. 
We find the equations and constants of motion through the Lagrangian formulation in the 
symmetric gauge. Our results are conveniently presented in dimensionless variables. In the 
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present view, bound states are obtained for a periodic potential with initial conditions, which 
satisfy Eq. lfTBj) and Eq. lfTTj) . We are able to write a simple expression for the frequency of 
motion, LUp = in the approximation of small oscillations. Trapped motion is obtained 

as a particular case of the motion of the center of mass with an additional restriction, which 
satisfy Eq. (j29|) . where there are no turning points. The definition of trapped states leads to 
a special range of values of the constant of motion, this is |T| < 1. 
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